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Abstract. 

Aims. We have investigated the accuracy and reliability of six methods used to determine the length of stellar bars in galaxies 
or N-body simulations. 

Methods. All these methods use ellipse fitting and Fourier decomposition of the surface brightness. We have applied them to 
N-body simulations that include stars, gas, star formation, and feedback. Stellar particles were photometrically calibrated to 
make B and K-band mock images. Dust absorption is also included. We discuss the advantages and drawbacks of each method, 
the effects of projection and resolution, as well as the uncertainties introduced by the presence of dust. 

Results. The use of N-body simulations allows us to compare the location of Ultra Harmonic Resonance (UHR or 4/1) and 
corotation (CR) with measured bar lengths. We show that the minimum of ellipticity located just outside the bulk of the bar is 
correlated with the corotation, whereas the location of the UHR can be approximated using the phase of the fitted ellipses or 
the phase of the m = 2 Fourier development of the surface brightness. We give evidence that the classification of slow/fast bars, 
based on the ratio R = RcR/Rbar could increase from 1 (fast bar) to 1.4 (slow bar) just by a change of method. We thus conclude 
that one has to select the right bar-length estimator depending on the application, since these various estimators do not define 
the same physical area. 

Key words. Galaxies: evolution - Galaxies: kinematics and dynamics - Galaxies: spiral - Galaxies: structure - Methods: 
N-body simulations 


1. Introduction 

Bars are ubiquitous and their importance for the long¬ 
term evolution of galactic di scs is well established nowa¬ 
days (e.g. see reviews by ISellwood & Wilkinson! Il993[ 
iKormendv & KennicutdE()o3 and references therein). One of 
the long-standing problems remaining for stellar bars is how to 
observationally delimit the area of a bar, i.e. to find where a bar 
morphologically ends or, in other words, how to measure its 
shape, length, and/or its width. Indeed, determining many cru¬ 
cial observational quantities relies on the definition of the bar 
size. However, a bar is not just an add-on morphological struc¬ 
ture embedded into an axisymmetric disc-Hbulge background. 
In fac t, some bulges (ornseudo-buIges. lKormendv & Kennicutd 
I 2 OO 4 I) could be the result of the secular evolution of bars, 
whereas part of the bar’s stellar material spends a significant 
fract ion of time far out into the disc (the so-called ’hot’ po pula- 
tion. lSnarke & Sellwoodll98^IPfenniger & Friedlill99ll) . This 
is the main reason any attempt to decompose ID photometric 
profiles or 2D surface brightness in a number of elementary 
contributions (bulge(s), disc(s), etc.) to recover the ’real’ bar 
by substraction, in most cases, is doomed to failure. 
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However, one of the simplest quantities used to estimate the 
strength of a bar is its ellipticity. For instance, lMartinl(ll99.fih 
found a correlation between the bar-axis ratio measured on 
blue photographic plates a nd star-formation activ ity in a sam¬ 
ple of 136 barred galaxies. Ichanelon et alJ lll999ll confirm this 
correlation using a more reliable photomet ry on red CCD im¬ 
ages. U sing a sub-sample of Martin’s data. iMartinet & Friedlil 
(I 1997 I) emphasize that the most actiye galaxies haye the longer 
and thinner bars, but they also conclude that this is only a nec¬ 
essary condition, not a sufficient one. Howeyer iKnaoen et al.l 
(l2002l) argue that this result could be dependent on the tech¬ 
niques used to determine bar lengths and strengths. Moreoyer, 
the bar strength esti mator Oh depends both on the axis ratio and 
the mass of the bar dSanders & Tubbsl[l98^ . 

There is another set of studies that needs accurate bar- 
length measurements. Numerous authors haye tried to corre¬ 
late the size of morphological structures to dynamical reso¬ 
nance locations. Circumnuclear and outer rings seem to be cor¬ 
related with the location of, respectiyely, the inner Lindblad 
resonance (ILR) and th e outer Lindblad resonance (OLR) (cf. 
iButa & CombesI [79961) . The ratio of the nuclear bar length 
to the large-scale bar length could be similar to the ILR-to- 
corotation ratio, the nuclear bar corotation being dynamically 
coupled with the large-scale bar ILR (e.g. iRaiitiainen & Salol 
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Il999l) . Howev er, some other sim ulations did not show such 
coupling (e.g. iHeller et alJEopl . In short, this matter is still 
being debated. 

The extent of the bar (either nuclear or large-scale) should 
thus be compared to the c orotation lo cation. There are sev¬ 
eral observational evidences (lKenlll990l and references therein) 
that stellar bars must end before the corotation. On the the- 


oretical side, the theory 
drodynamic simulations 

of orbits dContonoulosl 198f)l). hv- 

Sanders & Tuhh4l98f)l Athanassoulal 

Il992bl 

Regan & Teuben 

2004 etc.), and N-body simulations 


the corotation radius to the bar length is 1.2 + 0.2. Physically, 
this limitation is due to an increase in the amount of chaotic 
orbits close to the corotation. Thus, near the corotation, there is 
no orbital support to prolongate the shape of the bar. However, 
it is difficult to observationally conhrm this value since it 
not only depends on the method used to determine the pat¬ 
tern speed, but also on the criterion used to determine the 
bar length. For instance, lAguerri et al.l (l2003ll used four dif¬ 
ferent criteria to determine the bar length. Recent studies of 
the interaction between a stellar disc and a live dark halo have 
widely discussed the evolution of the ratio of the corotation ra¬ 
dius over the bar length, implicitly outl ining the importance of 
an ac c urate bar-length determination dOebattista & SellwoodI 
2222 L A^h^n^ss^ulaj^Misiri^iy. 200 j lAthanassoiilal l2003[ 


la & Misiri __ _ 

lO’Neill & Duhinskil2003llValenzuela & Klvpinl2003l) . 

Another difficulty arises from the morphological difference 
between mass and multi-wavelength surface-brightness distri¬ 
butions. Indeed, with N-body simulations, the bar length is de¬ 
termined on the mass distribution, whereas observational deter¬ 
minations use wavelength dependent surface-brightness distri¬ 
butions. Comparisons between the two imply an understanding 
and an estimation of the no-linear bias between mass and light. 
It poses, among others, the question of the reliability of the 
various criteria used in the literature and what they physically 
determine, depending on the support (mass or light) on which 
they are applied. 

To test the efficiency, reliability, and accuracy of vari¬ 
ous criteria for determining the length of a bar, we decided 
to apply them to N-body simulations that include stars, gas, 
and star-formation recipes. These simulations were photo- 
metrically calibrated, taking dust extinction into account (cf. 
Michel-Dansac & Wozniaki2004l hereafter Paper I). Thus, we 


could perform a systematic comparison of methods applied 
to mass distributions, as well as to photometrically calibrated 
mock images in blue and near-infrared wavelengths. Then, the 
location of the bar ends, estimated using six different tech¬ 
niques, was compared to the location of dynamical resonances. 
We thus exploit the comprehensive knowledge of the dynami¬ 
cal properties of the simulation to hnd the best bar-length esti¬ 
mator. 

In Sect. IJlwe describe the six criteria used to measure bar 
lengths. We also discuss some of their relative advantages and 
drawbacks. Then, in Sect. 0 we describe our numerical mod¬ 
els and their evolution. We also briefly recall the technique of 
photometric calibration used in Paper I. The six criteria are ap¬ 
plied to our simulations in Sect. 0] We also analyse inclination 
effects on our results in this section. Sect.|3deals with the tem¬ 


poral evolution of bar lengths. Correlations between the various 
bar length determinations and the location of dynamical reso¬ 
nances are studied in Sect. | 6 l The distinction between fast and 
slow bars and the limitations due to the lack of dark halo in 
our simulations are discussed in the same section. We Anally 
summarise our results in the last section. 

2. Measurements of bar length 

The method adopted bv lMartiiilll 19951) for measuring the length 
and width of bars is exclusively visual and relies on photo¬ 
graphic prints in the blue band. He estimates that the uncer¬ 
tainty is about 20% with this method. He defines the half-major 
axis of the bar as the length from the galaxy centre to the sharp 
outer tip where spiral arms begin, and the half-minor axis as 
the length from the centre to the edge of the bar, oval, lens, 
or spheroidal component, measur ed perpendicula r ly to th e ma¬ 
jor axis. Using CCD red images. Ichanelon et alJ ill999ll made 
a more sophisticated numerical analysis. They first extracted 
a 3-pixel wide photometric profile along the major axis of the 
bar. The half-major axis of the bar is the distance from the cen¬ 
tre of the galaxy to where the bar obviously ends. This is where 
the surface-brightness profile changes slope abruptly to become 
steeper; this also coincides with the origin of the spiral arms. 
This measurement is not automatic, as it relies on a subjective 
judgement of where the bar ends, and is comparable to that of 
Martin (1995). 

Beyond eye estimates, several attempts to And an automatic 
and objective criterion were made to determine where the bar 
ends and the disc or the spiral arms start. The objective criteria 
can be classified in two groups defined by the tools that used; 
3 based on ellipse fitting (noted Ei to E 3 ), and 3 based on the 
Fourier analysis (noted Fi to F 3 ). We applied all six criteria on 
the mass distribution and on the calibrated images. 

From now on, we use the expression ‘bar radius’ instead 
of ‘bar length’ since the various criteria used in this paper deal 
with either half-major axis length or radius. Each measurement 
must obviously be doubled whenever the real ‘length’ or the 
‘diameter’ should be determined. 

2.1. Ellipses fitting 

In the past, the analysis of isophotal shapes using a func- 
tional form was done using either classica l ellipses (e.g. 
IWozniak & Pierc3ll99l[ |Wozniak et al.lll99^ or generalised 
ellipses llAthanassoula et^ ilrnii), which permits the descrip¬ 
tion of boxy-like isophotes. 

By fitting el lipses to the isophote s of a large sample of 
barred galaxies, IWozniak et alJ ll 19951) determined some gen¬ 
eral rules about the radial behaviour of both the ellipticity (de¬ 
fined hy e — 1 - b/a where a and b are the half-major and 
half-minor axis lengths respectively) and the position-angle 
(PA) of the isophotes. The ellipticity is minimum or vanishes 
at the centre, because of either seeing effects or a luminous 
spherical bulge; then it increases to reach a maximum of¬ 
ten at about the middle of the bar, and then progressively de¬ 
creases towards at the place where the isophotes should 
become axisymmetric (disc) in the face-on case. Of course, 
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is determined by the galaxy inclination and/or the real non- 
axisymmetric shape of the disc. The PA is constant along the 
bar and, in general, sharply takes another constant value for the 
disc isophotes. If a spherical bulge is present, e will only start 
to increase outside the bulge. 

Three criteria can be dehned using the results of ellipse 
htting. For the first criterion (Ei), the bar radius is the ra¬ 
dius where the ellip ticitv prohle reaches a ma ximum. It has 
been introduced by IWozniak & Piercj ill99lll . who studied 
a sample of SBO galaxies in optical wavebands. Although 
this criterion is very usefu l for automated measurements (e.g. 
Re gan & Elmegre^ll997l Jungwiert et alJ [19971 iLaine et ^ 


2 QO 2 I 1 because the maximum ellipticity is always clearly de¬ 


hned in SBO galaxies, it seems to give rather short values of 
bar radius when it i s comp ared to eye estimates (cf. also Fig.0. 

IWozniak et alJ ( 1 19951) later introduced another criterion 
(hereafter called E 3 ) given by the position where the PA vari¬ 
ation exceeds 5° and ellipticity drops to disc values (rounder 
isophotes). This criterion gives higher valu es that were ex¬ 
pected to be upper limits of real bar radii. lErwin & Sparkel 
1I2OO3I) used a PA variation of 10° instead of 5, which should 
lead to slightly longer bars. The regions where the E 3 criterion 
is satished are generally less dusty than for Ei since gas den¬ 
sity strongly decreases outside the bar. The bar-radius estima¬ 
tion using E 3 might be more robust for most galaxies than Ei 
with respect to extinction but, depending on the detailed spatial 
distribution of the gas, could also lead t o large errors. A va ri- 
ant of the E 3 criterion has been used by I.Togee et HI (l2004l) to 
analyse almost 1 500 barred galaxies in the ACS GEMS survey, 
showing its reliability in data-mining studies. 

We now dehne a new criterion (E 2 ) that is also based on 
ellipse htting. In the region of the bar, where the ellipticity pro¬ 
hle looks like a bump, the PA is approximately constant, hence 
the PA prohle forms a plateau (cf. Pig. [I). Por the E 2 criterion, 
bar radius is measured at the end of this plateau. This corre¬ 
sponds to a twist between isophotes of the bar and isophotes of 
the disc region where spiral arms begin. This criterion gives an 
estimation of bar radii between those of Ei and E 3 criteria in 
the absence of dust extinction. 

To illustrate how these criteria could lead to different bar- 
radius estimations, we made an artihcial image that contains 
a bulge, two embedded bars, and a disc. Our purpose is not 
to make an extensive study, varying each parameter to quan¬ 
tify the quality of each bar-radius estimator. The bulge has a de 
Vaucouleurs prohle of central intensity - 9000, scale length 
Re - 200. It is slightly elongated along PA = 0° with an el¬ 
lipticity of 0.05. The exponential disc has a central intensity 
Id - 8000, scale length Rd - 300, and an ellipticity of 0.1 . The 
PA of the disc’s main axis is PA = 60°. The two bars have a 
surface density: 


l{x,y) = // 




3 ^ 

bi 


where a,- is the half-major axis (i.e. bar radius), h; the half- 
min or axis, /, the central inte nsity, and c, the shape parame¬ 
ter jAthanassoula et all Il990l) . The main bar has ai - 300, 
bi - 200, Ii = 5000, and ci varies from 2 at the centre (per¬ 
fect ellipses) to 3.5 at the ends, thus giving a rectangular-like 



radius (pixels) 


Fig. 1. Illustrative example of the prohles obtained by ellipse 
htting. Prom top to bottom: surface brightness in arbitrary 
units, ellipticity, and PA (in degrees) prohles as a function of 
half-major axis of the htted ellipse. Each point represents a 
htted ellipse. The vertical lines represent bar radii determined 
with El (dashed line), E 2 (dot-dashed line), and E 3 (dotted line). 


bar. It is aligned with the x-axis. A nuclear bar is added at 
PA = -45° with the same intensity prohle with parameters 
a 2 - 100, b 2 - 50, C 2 - 2 everywhere, and I 2 = 16000. All 
scale units are pixels and intensity units are numbers that could 
be considered as ADU. 

We display the results of the ellipse htting in Pig. [2 i.e. 
the ellipticity and PA of the isophotes as a function of the half¬ 
major axis length of the htted ellipse. Applying the above cri¬ 
teria leads to a bar radius of 255 for Ei, 280 for E 2 , and 310 for 
E 3 . This last criterion slightly overestimates the real bar radius 
because of the rectangular shape of the isophotes near the end 
of the bar. Ellipses are indeed not suited to strong rectangular 
bars. 
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Fig. 2. Illustrative example of the profiles obtained by Fourier 
analysis: bar-interbar contrast (top) and phase of the m - 2 
component (bottom) as a function of radius. Each point repre¬ 
sents a pixel. The dashed vertical lines symbolise the end of the 
bar determined using F 2 (top) and F 3 (bottom) criteria. 


2.2. Fourier analysis 

Fourier analysis of the surface brightness jOhta et al.[ Il99(]ll 
is based on the decomposition of azimuthal density profiles, 
I{r, 6), into a Fourier series 

00 

I(r, 6) = lo(r) + ^ [A,„(r) cos(m 6 ') -I- B„(r) sm(m0)] 

m=l 


where, for m 0 , the coefficients are given by: 
1 

A„, - — ( I(r, 6) cos(m6) dO 

^ Jo 
1 

Bm - — I lin ff) sin(m 0 ) dO. 

^ Jo 


The Fourier amplitude of the mth component (m > 0) is defined 
as: 

/„(r) = ^Al(r) + Bl(r). 

The bar-interbar contrast is then computed as: 


C - 4/fib - (4 + 4 + 4 + h)l(h - h + h- 4)- (1) 

We display C in Fig. |2l for the test image used above. 
Typically, C increases steeply and reaches its peak value in 


the bar regio n, then falls towa rd the bar end. The bar region 
is defined bv lOhta et al.l lll990l) as the zone where the contrast 
exceeds 2 (Fr criterion hereafte r). This criterion has been re¬ 
visited by lAgiierri et al.l il200f)l) . who redefine the end of the 
bar region where the contrast reaches the full width at half 
maximum (criterion F 2 ). In other words, the bar ends where 
C - 0.5 [max(C) + min(C)]. 

Another criterion (F 3 ) is determined using the radial profile 
of the phase of the m-2 component, which is defined as: 


(p2{r) = arctan[B2(4/A2(r)]. 


We display this phase profile in Fig.|2lfor the test image. The 
m-2 phase can be interpreted as the PA of the features 
that dominate the m-2 component. Thus, the phase suc¬ 
cessively gives the PA of the bar, then of the two-arm spiral 
structure, if present, or/and the PA of the disc (if not perfectly 
disky, which is always the case, due to inclination at least). The 
m - 2 phase profile is then similar to the PA profile, and the 
phase in the bar region is approximatively constant. We mea¬ 
sure bar radius at the end of this plateau. This criterion has 
already b een applied to mass distributio n from N-body simula¬ 
tions (e.g. lPebattista & Sellwoodl2000l) . as well as to observa¬ 
tions (e.g. lAgiierri et alJ2003h . 

These criteria applied to the artificial image lead to a bar 
radius of a; 280 for F 2 , whereas F 3 gives a; 270. However, Fi 
fails to give any estimation, since the contrast never exceeds 
the value of 2. Thus, for the case of our test image, F 2 and F 3 
underestimate the radius of the main bar. 


3. Numerical models 

3.1. hydro+N-body simulations 

We used PMSPHSF, the N-body code developed in Geneva. 
It includes stars, gas, and recipes to simulate star formation. 
The broad outlines of the code are as follows: the gravita¬ 
tional forces are computed with a particle-m esh method using 
a 3D polar grid jPfenniger & Friedlil Il993l) . the hydrodynam¬ 
ics equations are solved u sing the smooth particl e hydrodynam¬ 
ics (SPH) tec hnique fseelFriedli & Benzlll993l for this imple¬ 
mentation and lMonaghanll992l for a review of the method), and 
the star-formation process is based on Toomre’s criterion fo r 
the radial instability of gaseous discs jFriedli & Benzl 119951) . 
For the present work, the radiative cooling of the gas was com¬ 
puted assuming a solar metallicity for Runs A and B and a cos¬ 
mological metallicity for Run C. The metallicity chan ge of gas 
particl es was computed assuming net yields given bv lMaed^ 
At birth, star particles are created from gas, so they 
have the same metallicity. 

An initial stellar population was set up to reproduce a disc 
galaxy with an already formed bulge. These particles form 
what we call hereafter the ‘initial population’, as opposed to 
particles created during the evolution (‘new population’). The 
initial stellar positions and velocities for Ns particles of the 
same mass are drawn from a superposition of two axisymmet- 
ric Miyamoto-Nagai discs of mass Mi and M 2 , of scale lengths 
ai -i- b and a 2 H- b, respectively, and identical scale height b (cf. 
Table [Q where masses are in Mq and lengths in kpc). Their 
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Fig. 3. Evolution of the star-formation rate (SFR) of Runs A 
(black), B (red), and C (green). 

gravit ational potential is then given by jMivamoto & Nag^ 
Il975l) : 


where R — + y^. The gaseous component is represented 

by Ng particles for a total mass Mg, scale length Ug + bg, and 
scale height bg. As in Paper 1, we did not include a dark halo, 
to be able in the near future to compare the effects of having 
added such an additional component. 

Three simulations were used. Run A (cf. Fig. 0} was com¬ 
puted over 1 Gyr. It was analysed in some detail in Paper I. 
Runs B (Fig. 0 and C (Fig. were performed on a longer 
timescale. They will be used to check our results on the long 
term with different initial conditions. The main differences be¬ 
tween the three simulations are the initial mass of gas and its 
spatial distribution. The resulting star-formation histories are 
quite different (cf. Fig. |5}: Run A is more representative of a 
short timescale burst of star formation; Run C shows several 
small events but star formation is sustained over a longer pe¬ 
riod because of the more massive gaseous reservoir initially 
distributed over a greater length scale than for Run A; Run B is 
an intermediate case. 

During the three simulations, a bar progressively grows and 
then evolves secularly. The gravitational torques due to the 
bar and spiral structures drive the gas inwards. Star formation 
occurs along the transient spiral structure, along and around 
the bar, and in the central parts. The bar formation is thus 
associated to the starburst seen in Fig. 0 The bar-formation 
timescale differs from one simulation to the next, as a conse¬ 
quence of the different initial conditions. Another consequence 
is the morphological difference in the resulting bars. With these 


simulations, we have thus obtained a collection or ’sample’ of 
snapshots with a wide variety of bar shape, as can be seen in 
Figs. 1313 and 13 note the larger field of view used in Fig. |3 
than in the other two. Moreover, our simulated bars have prop¬ 
erties (such as ellipticity or strength) representative of nearby 
galaxies, as discussed in Sect.0] 

3.2. Mock calibrated images 

Dynamical properties of a bar obviously depend on the mass 
distribution. Collisionless N-body simulations are thus one of 
the best tools for studying them. However, when one deals with 
determining observational properties (i.e. photometric prop¬ 
erties in our case), an additional dependency makes the use 
of collisionless N-body simulations much more difficult; the 
mass-to-light ratio (M/F). Indeed, M/F varies in time and 
space since it depends both on the star-formation history and 
on the mixing of stellar populations by secular evolution. 
Simulations with gas and star formation permit photometri¬ 
cally calibrated light distributions to be obtained by performing 
a post-processing calibration with stellar-population synthesis 
models, as described hereafter. 

For each stellar particle, given its age and metallicity, the 
mass-to-light ratio was obtain ed from a bi-linear interpo lation 
of the tables of GISSEF2000 dBruzual & Chariot! 1 1993I) . for a 
Salpeter initial mass function with a lower mass cutoff Mj - 
0.1 Mq and a upper cutoff = 100 Mq. 

For the initial population of Run A (present at the begin¬ 
ning of the simulation), we assumed an age of 10.4 Gyr and 
metallicity of Z = 0.004, like the model A of Paper I. In that 
paper, we considered various calibrations of the initial popula¬ 
tion and found that the length of the bar (determined with the 
E 3 criterion, defined above) is almost independent of the cali¬ 
bration model. We thus decided here to restrict our study to one 
calibration model for the sake of clarity. 

Obviously, the assumptions on the age and metallicity of 
the initial population implicitly imply that all particles were 
simultaneously born 10.4 Gyr before the beginning of the sim¬ 
ulation. Assuming an observer located at z = 0 at the end of 
the simulation (f = 1 Gyr), this calibration implies a redshift 
of formation of z = 3. For the longer Runs B and C, we consid¬ 
ered several sets of initial population ages and metallicities. For 
clarity, we deal only with the calibration model that assumes an 
initial population age of 4.4 Gyr and a metallicity of Z = 10 
This leads to a final age of 11.4 Gyr, the same as for Run A. 

In Paper I we worked with four wavebands to cover the 
wavelength domain from the visible to the near infrared, 
namely B, R, H, and K Johnson. Here, we restrict our analy¬ 
sis to B and K bands without any loss in generality. 

Figures |3 to |3 display the evolution of the mass density 
and the B and K calibrated images with and without dust, in 
the bar region for a selection of snapshots. For Run A, large 
scale snapshots of the evolution can be found in Paper I. These 
figures clearly show that the disc undergoes a complete redis¬ 
tribution of the mass. For instance, for Run A after 1 Gyr, the 
total mass inside the central kpc has roughly doubled. The first 
cause of this mass inflow is the formation of a stellar bar. Then, 
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Table 1. Initial parameters of the simulations 


Run 


N, 

Ml 

M2 

M, 

Q,[ 

Cl2 

b 

Ug 

bs 

A 

500000 

50000 

"1(F 

10" 

1.1 10"’ 

0.5 

3.0 

0.5 

3.0 

1tP~ 

B 

500000 

10000 

lO*" 

10" 

3.66 10*“ 

0.5 

3.0 

0.5 

3.0 

0.5 

C 

600000 

10000 

210*“ 

10" 

410*“ 

0.5 

5.0 

1.0 

10.0 

0.5 









Fig. 4. Evolution of Run A from f = 0 to t = 1 Gyr. From left to right are displayed the mass distribution (in log of Mq pc 
and the calibrated images (B dust-free, B with extinction, K dust-free, and K with extinction) in mag arcsec^^. The field of view 
(10 kpc) is the same for each frame. The particles have been rotated so that the bar position is roughly horizontal. Greyscale 
images and isocontours of mass surface density range from lO' to 10^ ^ Mq pc“^. Mass isocontours are logarithmically spaced 
by 0.47 log(MQ pc“^). The B greyscale images range from 11 to 26 mag arcsec^^, while contours range from 20 to 23 mag 
arcsec^^ and are spaced by 1 mag. The K images range from 9 to 22 mag arcsec“^; K isocontours range from 17 to 20 mag 
arcsec^^ and are spaced by 1 mag. 



Michel-Dansac & Wozniak: The length of stellar bars 


7 


Fig. 5. Like Fig. 0]but for Run B from f = 0 to f = 5 Gyr. Greyscale images and isocontours of mass surface density range from 
10' to 10^ ^ Mq pc“^. Mass isocontours are logarithmically spaced by 0.47 log(MQ pc“^). The B greyscale images range from 
12 to 27 mag arcsec^^, while contours range from 19 to 23 mag arcsec^^ and are spaced by 1 mag. The K images range from 10 
to 23 mag arcsec^^; K isocontours range from 16 to 20 mag arcsec^^ and are spaced by 1 mag. 


Fig. 6 . Like Fig. EJbut for Run C from f = 0 to f = 5 Gyr. The field of view is 16 kpc. Greyscale images and isocontours of 
mass surface density range from 10' to 10^ Mq pc“^. Mass isocontours are logarithmically spaced by 0.44 log(MQ pc“^). The B 
greyscale images range from 13 to 27 mag arcsec^^, while contours range from 19 to 23 mag arcsec^^ and are spaced by 1 mag. 
The K images range from 10 to 23 mag arcsec^^; K isocontours range from 16 to 20 mag arcsec^^ and are spaced by 1 mag. 


due to the gravitational torques exerted on the gas by the stellar 
bar, the extra mass in the form of gas and new stars amounts 
to 3.5 10^ Mq at f = 1 Gyr for Run A, which is only 30% of 
the whole additional mass. In fact, the initial stellar population 
contributes to the other 70%. 

After individual particle photometric calibration, mock 
CCD images were obtained summing particle luminosities into 
a 512x512 pixels grid. The field of view was 60 kpc, which 
gave a spatial resolution of 117 pc, 1.3 times our small¬ 
est A-body grid resolution. We thus produced one frame per 
waveband and per snapshot of the simulation. Our results are 
obviously independent of the bar PA with respect to the North 
or any other axis. We thus decided to systematically rotate the 
positions of particles to align the bar with the x-axis. 

To mimic real observations we should have to convolve our 
images with a point-spread function. However, this last stage 
depends on the telescope and observation-site characteristics. 
It thus introduces a few free parameters that cannot be con¬ 
strained without any detailed comparisons with real observa¬ 
tions, which is not our purpose. 

3.3. Dust extinction 


The ellipticity profiles in both bands are significantly dif¬ 
ferent for f < 1 Gyr. In Paper I, we proposed that, even in the 
absence of extinction, two causes are simultaneously responsi¬ 
ble for the wavelength dependence of ellipticity profiles. First, 
the spatial distribution of the new population is very elongated 
along the bar because of new stars born in the gas flow along 
the bar that is narrower than the stellar bar. Then, secular evo¬ 
lution is responsible for progressively making it rounder. This 
explains the high ellipticity reached in the bar when the star- 
formation rate (SFR) is high (f < 0.5 Gyr for Run A) and, 
only in part, the subsequent decrease. Second, the luminosity 
ratio between the new and the initial populations is wavelength 
dependent, being higher in B than in K band. After the SFR 
maximum, the morphology becomes gradually dominated by 
the luminosity of the initial population that has a rounder spa¬ 
tial distribution. When dust extinction is taken into account, 
especially in the B band, there is no longer a unique be¬ 
cause the real maximum is located in the dustiest region (e.g. 
f = 0.3 Gyr in Fig.|3i. A comparison between B and K-band 
El measurements confirms that ellipticity is strongly dependent 
on the colour. 

It is interesting to compar e the maximum ellipticity with 
the bar strength Qh defined as jCombes & Sande^ll98lh : 


Dust extinction in B and K bands was simulated assuming a 
constant gas-to-dust ratio 

A(Hi)/Av = 5.34 X lO^'cm'l 

Then, Ay is converted to A^ and A^. Extinction was computed 
in a cube with the same spatial resolution as for our images 
and 11 slabs along the line-of-sight (cf. Paper I for more de¬ 
tails). Each slab absorbs the stellar luminosity behind it. Eor 
each slab, the gas density distribution is obtained by convolv¬ 
ing particle positions with the SPH kernel. 


4. Tests of bar-radii criteria on numericai 
simuiations 

Each mock image was analysed by fitting simple ellipses to 
the isophotes. Profiles of surface brightness, ellipticities, and 
PA were obtained by increasing the half-major axis length a 
by a factor 1.01 between each fit. We chose such a low value 
to obtain a good resolution in the inner region. We display the 
results for Run A in Pigs. 0 and| 3 for the same selected times as 
in Pig.0 Results for the other runs are qualitatively similar, so 
they are not displayed. 


Qb = max 


77max 
^0 


(r) 


< Fair) > 


where is the maximum tangential force at radius r and 

< > is the average radial force from the axisymmetric 

component. Recently, this bar strength estimator has been also 
used in photometric studies ilButa & Blockl 1200 ll iBiita et all 
l2005t) . Comparison of the evolution of the bar ellipticity with 
that of Qb (cf. Pig. [TOJ clearly confirms that the bar axis ra¬ 
tio (or ellipticity) is not an estimator of the bar strength Qb. 
Indeed, bars with larger quadrupole moments and lower axis 
ratios can have the same Qb as bars with smaller quadrupole 
moments and higher axis ratios. Moreover, the radius at which 
Q(r) reaches its maximum, used to define Qb, is always located 
in the circumnuclear region. Thus, Q(r) or Qb is not useful for 
determining the bar radius. Th e distributions of O h for our sim¬ 
ulations and for real galaxies jRiita et al .1 l20n.5li are very sim¬ 
ilar, except for Qb < 0.15 (very weak bars). Bar morphology 
and properties, such as strength and ellipticity, produced in our 
simulations are thus representative of real galaxies in the local 
universe. 

Due to the definition of E 2 and E 3 , PA variations are 10° 
within the radial range defined by E 2 and E 3 criteria. This region 
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Fig. 7. In each panel, for the two bands (B plotted in blue and K in red) and the mass distribution (in black) for Run A: from top 
to bottom, surface brightness (ju), ellipticity (e=l-b/a), PA, and central shift for each fitted ellipse. The range in radius for the 
figures begins at 1 px, i.e. 117 pc; the upper bound corresponds to a region that encloses the bar during the entire run. The mass 
density profiles have been scaled so that the 30'^ magnitude corresponds to 1 Mq pc“^. The labels Ei, E 2 , and E 3 refer to the 
criteria used to determine the bar radius in B band. See the text for their definitions. 


of strong isophotal twist clearly marks the boundary of the bulk 
of the bar. Outside this transition region the surface photometry 
is dominated by the disc contribution and possibly by spiral 
arms; however, the size of this twisted region could be very 
time dependent, especially in the presence of spiral arms. This 
leads to variation from one snapshot to the next in E 2 and E 3 
bar radii since the region could shrink or expand on a short 
timescale. 

Fourier analysis was also done for each mock image. Radial 
profiles of the contrast computed with B and K images, as 
well as with the mass distribution, are displayed in Fig. for 
Run A. In general, the contrast follows a well-defined curve, 
with a clear maximum once a bar is formed. However, dur¬ 
ing the most active star-formation phase, large fluctuations of 
the contrast, especially in the dust-free B band, make it diffi¬ 
cult to draw any radial profile. The contrast is also very sensi¬ 


tive to the extinction in the B band, as opposed to the K band. 
Interestingly, the contrast is always higher in the B band than in 
the K band. The lowest values are reached for the mass distri¬ 
bution. With another observational tool, this confirms the result 
found in Paper I that the bar shape (ellipticity or contrast) de¬ 
pends both on the wavelength and on the delay after the major 
star-formation phase. 

However, it could be difficult or even impossible to define a 
maximum (e.g. t - 0.3 or f = 0.4 Gyr in Fig.|3. Moreover, the 
minimum could sometimes be negative in the case of very elon¬ 
gated and thin-light distribution, which is the case when star 
formation is active along the gas streamlines in the bar. Indeed, 
/ib takes negative values. This also explains why we failed to 
apply the Fi criterion to dust-free B images. Sometimes, the 
contrast remains lower than 2 everywhere (as for the test image 
in Sect. which makes it impossible to obtain F 1 . 
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Fig. 8 . Like Fig. □but for dusty images. Results for the mass density are the same than in Fig. □since they are not affected by 
extinction. 


4.1. Projection and resolution effects 

Because real galaxies are never exactly face-on, it can be asked 
how the bar radius determined is affected by projection effects 
when using the six criteria. We used one snapshot (Run A at f = 
0.77 Gyr) to study, on one hand, the effects of spatial resolution 
and, on the other, the effects of inclination. We used images in 
K band. 

When working with inclined galaxies, we also had to take 
the PA of the bar with respect to the line of node into account. 
Therefore we measured the radius of the bar using the six crite¬ 
ria for simulated galaxy images inclined from 15° (nearly face- 
on) to 70°. For each inclination, we used 4 values for the bar 
position-angle (PA = 0°, 30°, 60°, and 90°). 

For i - 15°, the measurements of bar radii are rather insen¬ 
sitive to the PA variation. For i > 15° there are two regimes. 
For PA < 60°, bar radii decrease with PA, and this decrease 
is stronger for higher values of inclination. For PA > 60°, the 
trend depends on the criterion. For example, bar radii deter¬ 
mined with Fi and F 2 decrease with PA, whereas they increase 


for E 2 . With the Ei criterion, the bar radius is almost constant 
in this range of i and PA. For E 3 and F 3 , bar radii depend on the 
degree of inclination. In addition for i = 70° Fi does not work, 
and for high values of i and low values of PA, Ei gives a bar 
radius close to the values obtained with E 2 , F 2 , and F 3 . 

To give a rough estimate of the error made on the bar length 
using imperfectly deprojected images of inclined galaxies, we 
find that in all cases the difference for i < 30° is less than 20% 
and the difference for i < 15° is less than 10%. For our images, 
with a bar radius at 5 kpc, 10% represents less than 5 pixels. 

Finally one could also wonder how the spatial resolution 
of the images affects the determination of the bar length. We 
then measured the bar radius using those six criteria in K-band 
images with lower resolution (two to hve times lower reso¬ 
lution; the lowest resolution is then approximatively 0.6 kpc 
per pixel). For a two-times lower resolution, there is no dif¬ 
ference: variation lower than 2 %, which is lower than incerti¬ 
tudes due to human measurement. For lower resolution, bar- 
radius measurements with E 2 , Fi, and F 2 have small variations 
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T = 0.2 T = 0,3 T = 0.4 



Radius (kpc) 



Fig. 9. Contrast profiles inside the central 6 kpc of Run A, determined on the mass distributions (black crosses), B images with 
extinction (blue full circles), dust-free B images (blue line), and K images (red line). Measurements inside the central kpc are not 
plotted since they are not reliable. 


(=s 5%), whereas with Ei, E 3 , and F 3 they vary more signifi¬ 
cantly (10- 15%). When resolution is decreased, the trends are 
that bar radius determined with £3 decreases slightly, increases 
slightly with Fi and F 2 , increases with E 3 , and decreases with 
F 3 . The El criterion is more fluctuating, and bar radius is then 
smaller or longer, depending on the resolution. Due to these ef¬ 
fects, El, £ 2 , and F 3 criteria could give close values for the bar 
radius at low resolution. 


4 . 2 . Dust extinction effects 

Extinction effects in the bar region are important. For Run A at 
t » 0.3 Gyr, less than 20% of the bar B-band luminosity es¬ 
capes, increasing to 55% at f 1 Gyr. As shown in Paper I, 
B isophotes are more disturbed by dust extinction than K ones. 
Since our dust model is proportional to H i column densities, 
extinction decreases with time following the consumption of 
gas by star formation. A comparison of Figs. 0 and |3 shows 
that ellipticity profiles are completely different and very noisy 
when extinction is high, making bar-radius determination diffi¬ 


cult. Differences in B-band ellipticity profile are still present at 
f = 1 Gyr but mainly appear in the central kpc. 

Regions with high SFR are located where gas density is 
the highest. In these regions, extinction is then very high since 
it is proportional to the gas column density. The regions that 
contain the youngest and bluest populations are also the most 
obscured. Since the gas is consumed during the evolution and 
since gas and stars have different kinematics, the gas density 
decreases secularly, enabling the escape of the blue light af¬ 
ter several tens of Myr. However, the gas consumption occurs 
inhomogeneously, and at various rates, in the bar region. 

All criteria are not affected by extinction in the same way. 
For instance, the E 3 criterion is in general rather insensitive to 
the extinction because the measurements are made in regions 
where the gas density is lower than in the innermost regions 
but, depending on the detailed spatial distribution of the gas, it 
could also lead to large errors. For Run A, the E 3 bar radius at 
t - 0.6 Gyr differs by 1 kpc, whether extinction is present or 
not. It is also noteworthy that, although the Ei bar radii differ 
when extinction is added to the images, maximum ellipticity 
values do not significantly differ from the dust-free case. For 
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and E 2 in B band give similar results during the active star- 
formation phase (f ^ 0.6 Gyr) because of the steep ellipticity 
gradient due to the sharp bar boundary. 

The bar radii for the mass density and K images show the 
same trend as a function of time. They thus could be consid¬ 
ered as similar on average, although results clearly show larger 
fluctuations for the mass density than for the K images. 

For the three criteria, the large fluctuations of bar radius 
from one snapshot to the next do not have a dynamical origin. 
This is only the consequence of our attempt to apply the criteria 
in the most objective way. Indeed, we systematically disregard 
bar-radius determination for previous snapshots. These fluctu¬ 
ations must be considered as an estimation of the error caused 
by using a criterion based on ellipse fitting. Since an ellipticity 
maximum is easier to determine than any other peculiar place 
on ellipticity or PA profiles, Ei gives less noisy estimations than 
E 2 or E 3 . 


5.2. Results with Fi, F 2 , and F 3 


Fig. 10. Evolution of the bar strength Qt (open circles) for Runs 
A (top), B (middle), and C (bottom). The maximum ellipticity of 
the bar is plotted as full circles for the mass density (black), B 
images (blue points), and K images (red points). 

Run B, the ellipticity and PA profiles in B band are strongly 
affected by dust extinction during the 2 first Gyr. For most of 
these snapshots, we were not able to determine the bar radius 
using El, E2, and E3. 

Paradoxically, the contrast profiles C(r) are much more re¬ 
liable when extinction is taken into account, since the very ac¬ 
tive regions along the gas streamlines are hidden. Fluctuations 
of /ib, the denominator in the definition of C, are thus strongly 
smoothed. The application of F 3 to dusty images leads to much 
more difficult measurements of bar radii, in particular in the B 
band and during the first hundred Myr of evolution when the 
gas density is still high. 

5. Evolution of bar radius 

The evolution of the bar radius, determined using all six crite¬ 
ria, is plotted in Figs.^Jand^] In the case of Runs A and B, 
for t $ 0.2 Gyr, the bar progressively grows but is still too weak 
to give a sharp signature in the surface-brightness distribution. 
Thus, each of the criteria fails to accurately define a bar radius 
during this growing phase. This is also the case for Run C for 
f 1 Gyr, since the bar formation has a longer timescale due 
to the more massive initial bulge. 

We hereafter discuss the results obtained for Run A and then 
look into details of the evolution on the longer term. 

5.1. Results with Ei, E 2 , and E 3 

Criteria based on ellipse fitting give different estimations of 
the bar radius; Ei gives generally lower values, while E 3 gives 
the highest estimation, as expected (cf. Sect. O) However, Ej 


In comparison with the previous criteria, Fi, F 2 , and F 3 give 
different radius estimations, evolution, and level of fluctuations. 
Indeed, the bar radius determined with Fi and F 2 is roughly 
constant with a little dispersion for f > 0.5 Gyr, whereas it 
increases slightly with time and with a stronger dispersion for 
the other ellipse-based criteria. Bar radius determined with F 3 
first increases, then slightly decreases for f > 0.8 Gyr. For F 1 
and F 2 , the level of fluctuation from one snapshot to the next 
is lower than in the case of E 2 and E 3 because; 1 ) these criteria 
do not rely on human decision and 2 ) high frequencies in the 
image are smoothed out. However, F 3 is a bit noisier than Fi or 
F 2 probably because this criterion also relies on human decision 
as do the E 2 and E 3 criteria. Moreover, Fi and F 2 give similar 
results for the mass but differ on K-band images, F 1 giving a 
longer bar than F 2 by » 1 kpc. 

The F 2 criterion deserves particular attention since it gives 
similar and reliable results for the mass, B, and K images 
for t > 0.5 Gyr, that is, 0.2 Gyr after the peak of SFR. For 
t < 0.3 Gyr, when the bar is still growing, bar-radius estima¬ 
tion is unreliable. For 0.3 < t < 0.5 Gyr, when star formation 
is the most active, the bar radius is shorter in B than in K, the 
highest value being obtained for the mass. However, the con¬ 
trast C profiles (Fig. clearly show that these estimations are 
not reliable, apart from the mass density. 

Fet us recall that it was not possible to get valuable bar radii 
with F 1 on B-band images because the contrast often remains 
above the threshold (fixed at 2 by Ohta et al.) until the end of 
the bar. 

5.3. Long-term evolution 

All the above analysis were applied to the other two runs that 
were performed until t - 1 Gyr. We reached roughly similar 
conclusions when we compared the results of the three runs. 
Such a comparison does make sense only on condition that the 
different evolution histories of each run are taken into account. 
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Run A Run B Run C 



0.0 0.2 0.4 0.6 0.8 1.0 0 2 4 6 0 2 4 6 

Time (Gyr) Time (Gyr) Time (Gyr) 


Fig. 11. Bar and resonance radii for Runs A, B, and C. For each panel, lines represent the ILR (shorter radius), UHR (intermediate 
radius), and CR (greater radius). The dimmed region shows the Li, 2 -T 4,5 range. Bar radii determined with the maximum ellipticity 
criterion (Ei) are plotted as open triangles, those determined with the minimum ellipticity criterion (E 3 ) are plotted with open 
squares, and the new criterion (E 2 ) is represented by full circles. Top panels: bar radii determined for the mass. Middle panels: 
B dust-free images. Bottom panels: K dust-free images. 


The first Gyr of Run A is therefore comparable to almost the 
first 1.5 Gyr of Run B and the first 3 Gyr of Run C. 

The differences between Run A and the two other runs 
mainly concern the criteria Ei and Fi to F 3 . For Run B, the Fi 
criterion can be measured from B-band images, and the results 
are similar to those in K band. For Run C, the measurements us¬ 
ing this criterion are possible and reliable for f ^ 2 Gyr and give 
the same values as the F 2 criterion. Another difference concerns 


the El criterion for Run C. Indeed, it gives approximatively the 
same values as the E 2 criterion for f $ 3 Gyr. The F 3 criterion 
gives a significantly greater bar radius than those obtained with 
Fi and F 2 . This criterion is very comparable to E 2 . 

On the other hand, for the criteria E 2 , E 3 , and Fi to F 3 , each 
gives similar results when we determine bar radius on the B or 
K-band images or on mass distributions, taking the uncertain¬ 
ties of measurements into account. 
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Run A Run B Run C 



Time (Gyr) Time (Gyr) Time (Gyr) 


Fig. 12. Like Fig. 11 1 Ihiit for bar radii determined with the Ohta et al. criterion (Fi, open triangles), Aguerri et al. criterion (F 2 , full 
circles), and m - 2 phase criterion (F 3 , open squares). Top left panel: measurements for the mass. Middle panels: B dust-free 
images. Measurements for t < 0.4 Gyr are not reliable for Run A (see text). Bottom panels: K dust-free images. 


6. Discussion 

6.1. Comparison with dynamicai radii 


K\r) = 


520 520 

dx^ 5y2 



In Fig. CH we also plot the radius of the CR (Rcr), ILR, and 
UHR resonances. In order to estimate the dynamical reso¬ 
nances radii, we computed the cir cular orbit freque ncy Q. and 
the radial epicyclic frequency k as ilPfennigeilll990l) : 


Q 2 (r) = 


/ 1 ^\ 

\7'^/ 


where < 1 ) is the gravitational potential and (■ ■ ■) stands for an 
azimuthal average. 

Strictly speaking, these frequencies predict the oscillation 
frequencies of the orbits only in the axisymmetric case. They 
do not indicate whether families of periodic orbits actually fol¬ 
low such oscillations when the bar growth breaks the axisym- 
metry. In fact, a proper orbital analysis is more appropriate. 
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Fig. 13. Evolution of the relative difference between Li ,2 and 
L 4,5 (top) and bar strength Qb (bottom) for Runs A (full circles), 
B (lozenges), and C (squares). 

However, a number of previous orbital studies incline us to 
consider that the epicyclic approximation could lead to an ac¬ 
ceptable estimation of th e resonance locations. 

Indeed, for the ILRs, lAthanassoulal (Il992al) has computed 
the spatial extension of the X 2 and x^, families associated with 
the existence of these resonances. Her comparison with ILR 
radii calculated by averaging the frequencies over the az¬ 
imuthal angle shows that the error cannot exceed 10 % in the 
worst cases. Since our A-body mode l s are w ell inside the pa¬ 
rameter space used by lAthanassoulal (Il992al) . we confidently 
use the linear approximation to estimate the radius of the outer 
ILR. 

To our knowledge, the differences between the real 
Lagrangian points location and the axisymmetric corotation ra¬ 
dius have been less quantified in the past. In the case of barred 
galaxies, Lagrangian points are defined as the points of equi¬ 
librium between centrifugal and centripetal forces along the 
main axis of the bar perturbation. The radii of these points con¬ 
verge towards the c orotation c i rcle w hen the bar perturbation 
vanishes. Lollowing IPfeWigei] lll99oll and assuming a perfect 
bisymmetric barred potential, we call Li _2 and L 45 the radii 
of the Lagrangian points along the major-axis and intermedi¬ 
ate axis (minor-axis in a face-on projection), respectively, of 
the bar perturbation. In a strong bar potential, L 4_5 is slightly 
smaller than Li, 2 . However, the amplitude of the difference be¬ 
tween Li ,2 and L 4 5 should be roughly proportional to the bar 
strength. 

In Lig. [21 we thus display the quantity AL = 
(Li ,2 - L 4 , 5 )/L, where L is the average between Li _2 and 
The relative amplitude of the difference between L \^2 and L 4_5 
rarely exceeds 15% even in very strong bar phases, as for in¬ 
stance during the bar-formation phase in Run B. A standard 


value for a slowly evolving bar seems to be in the range 5 to 
10%. A comparison with the Qb evolution clearly confirms this 
expected trend. When the bar gets stronger, the difference be¬ 
tween L] 2 and L 4 5 increases. 

In the epicyclic approximation, corotation is the radius 
at which Q(Rcr) = where Q. is the circular orbit 

frequency and Q.b is the bar pattern speed. The value of 
Q.b was accurately determined using the Tremaine-Weinberg 
method ( iTremaine & Weinberg! 1 19^ since it has been shown 
that the method is very reliable on stell ar yelocity fields 
dHernandez, Wozniak, Carignan et alll2005l) . i.e. the same yal- 
ues than those obtained using the phase of the bar. In Lig. 11 Hwe 
compare the corotation radius with Li ,2 and As expected, 
the corotation radius lies between the Lagrangian points but is 
closer to Li _2 than Thus, we can use the linear approxima¬ 
tion to compute the corotation radius and AL, as defined aboye, 
could be used as an estimator of the error made on the corota¬ 
tion location. 

Keeping in mind a typical error of 10%, as discussed aboye, 
we haye determined the position of the radial ILR defined as 
the solution of Q(Rilr) - a-(Rilr)/2 = Q.b and UHR for which 
f^(^UHR) - M^UHr)/ 4 = Qi. 

The last check was to compute the orbital frequencies k and 
fl of a representat iye sample of particle s. We applied a yariant 
of the technique of lAthanassoulalll2003ll . We froze the potential 
at a giyen time in the simulation and then computed orbits in the 
inertial pote ntial to determine the main fr equencies using the 
technique of ICarointero & AguilaJ lll998l) . We used between 
100000 and 200000 particles as initial conditions chosen at 
random. This leads to the determination of k and Q. for each 
particle which can be plotted in a a- - Q diagram. Resonant 
families (e.g. ILR) and corotation are thus easily identifiable. 
We repeated this yery CPU time-consuming computation' for 
4 snapshots in order to check the consistency with the linear 
approximation described aboye. We thus determined that the 
error on the corotation radius neyer exceeds 10 %. 

The resonance locations are obyiously independent of our 
photometric calibration and thus could be compared to our es¬ 
timation of bar radii using yarious criteria. Howeyer, for all the 
reasons mentioned in Sect. |3 about Ei cannot be, a pri¬ 
ori, a good estimator of any dynamical radius (ILR, CR, etc.). 
This is clearly confirmed in Lig. CH Also, E 3 giyes bar radii 
that are yery close to the corotation. This strongly suggests that 
the ellipticity minimum at the end of the bar could be due to 
the signature of the corotation radius. Howeyer, if we define 
the corotation region as between Li _2 and E 3 giyes a lower 
limit of the corotation. 

The criteria E 2 , F 2 , and F 3 giye bar radii that are rather well 
correlated with the UHR radius. Howeyer, F 2 is less correlated 
when SLR is high (for t < 0.5 Gyr), whereas E 2 could giye a 
higher error on the UHR position. One more difference appears 
between yarious simulations for the F 2 criterion. Lor Run C, it 
is correlated with the UHR radius for t ^ 3.5 Gyr, i.e. for the 
duration comparable to Run A. On longer timescales, this is no 
longer the case. Lor Run B, the F 2 criterion neyer appears to be 

' It takes roughly two weeks by snapshot on a 12-node cluster for 
200 000 orbits 
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correlated with the UHR radius, and it gives values between the 
ILR and the UHR radius. Nevertheless, E 2 and F 3 criteria are 
still correlated with the UHR radius on longer timescales, like 
the E 3 criterion with the corotation radius. 

The fact that E2, F3, and in some cases F2 criteria are 
good estimators of the UHR radius during the hrst Gyr can 
be understood in the framework of s tellar orbits. Indeed, it 
has been shown jContonoulosl Il981all that the UHR (or 4/1 
resonance) is a gap where no family of periodic orbits can 
exist. Moreover, between the UHR and the corotation, most 
families of periodic orbits are unstable jContopoulosl Il981bh 
if the bar perturbation is stronger than 10 % of the potential 
backgr ound (the stellar disc). Using 3D models. Fskokos et'afl 
1120021 ) also conclude that the most appropriate orbits to sustain 
a bar are those inside the UHR. Thus, between the UHR and 
the corotation, a bar is mostly sustained by semi-chaotic orbits 
llWozniak & Pfennigerl Il999l) . The density response of these 
semi-chaotic orbits in the conhguration space is rounder than 
the response density of orbits trapped around the major fam¬ 
ilies of periodic orbits (i.e. xi and 3/1 resonant family). Thus 
the contrast as dehned by Eq. Q should strongly decrease af¬ 
ter the UHR. This also explains why E 3 is a good tracer of the 
corotation radius. 

However, for Run C, the F 2 criterion is no longer correlated 
with the UHR radius from f ^ 3.5 Gyr and never appears to 
be correlated with the UHR radius for Run B. This means that 
this criterion is not a reliable estimator of the resonance loca¬ 
tion. The bar-interbar contrast indeed shows the pattern of the 
bar for all the snapshots of these three runs, but the decrease 
of the contrast at the end of the bar, which is the transition re¬ 
gion between the bar and the disc, seems to be very sensitive 
to local conditions (such as star formation at the end of the bar, 
the spiral structure in the disc, and the thickness of the disc); 
therefore, these measurements depend on the type of bar and 
evolution. Then, even if the F 2 criterion is an automatic and re¬ 
liable criterion of the bar radius, it dehnes a radius that is not 
always correlated with a resonance and then is not comparable 
from one to another galaxy. 

6.2. Fast and slow bars 

We plotted the corotation versus the two more accurate deter¬ 
minations of bar radius (cf. Fig. HU). The corotation radius was 
estimated using L - (Li ,2 + T 4 , 5)/2 as in Sect. 16. II The error is 
given by the difference between Li 2 and L 4 , 5 . 

Using such a plot, lAgiierri et alJ ll2nn3h have divided fast 
bars (1 < 'R - RcR/Rbm < 1-4) from slow ones (R > 1.4). 
For our simulations, especially for the long Runs B and C, R 
does not evolve very much with time except during the bar- 
formation phase. Once the bar is settled in, R remains rather 
constant. 

It is noteworthy that the value of R is very sensitive to the 
dehnition of the bar radius. For the same simulation (either 
Run B or C), E 2 (or F3) leads to the conclusion that the bar 
is a rather slow rotator, while E 3 leads to the opposite conclu¬ 
sion. There is unfortunately no obvious theoretical dehnition 
of the ‘radius’ or ‘length’ of the bar. The only clean theoret- 



Fig. 14 . Plot of corotation radius versus bar radius for Runs B 
(lozenges) and C (squares). The corotation radius was estimated 
using L, and the error is given by the difference between Li _2 
and L 4 S. For the bar radius, we used criteria E 2 (full symbols) 
and E 3 (open symbols), determined on K-band images. 


ical predicti on is that bars cann ot extend outside the corota¬ 
tion radius dContopoulosl Il980l) . whereas they could be lim- 
ited by the disc scale length and thus could end near the IFR 
dCombes & Elmegreenlll993l) . Moreover, when the ratio of the 
bar radius to the coro tation radius is given (e.g. 1.2 + 0 . 2 , cf. 
lAthanassoulall 19925) . the bar radius is in fact dehned as the 
half-major radius of a Ferrers ellipsoid that is not an observ¬ 
able radius, so that any comparison with observational bar radii 
could be difhcult. 

In the case of the measurement of R, the use of E 3 should 
obviously be avoided. Indeed, E3 points out the corotation, 
leading to a ratio R ^ 1. Thus, E2 and F3 criteria, being loosely 
correlated to a dynamical radius (UHR), could be useful in a 
hrst estimation of a bar radius. Dehning the bar radius in such 
a way could be considered as a loose dehnition of the bar radius 
since the UHR does not limit the extension of the bar. We can 
thus imagine converting Ruhr into but the ratio Rcr/Ruhr 
decreases from 1.5-1 .8 during the short bar-formation phase to 
1.3 when the bar is settled in. Thus, the location of the coro¬ 
tation could only be retrieved within an error of the order of 
10 %. Of course, the combination of the E2 criterion with E3 
could lead to a better estimation of the corotation radius. 

6.3. Dark-halo effects 

The main effect of a live dark halo (except to hatten the ro¬ 
tation curve of the disc at large distance) is to permit ex¬ 
changes of angular momentum with the stellar disc. The rate 
and amount of these exchanges depend on the velocity disper- 
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sion of b oth the disc and the ha l o, and on the relative halo 
mass (cf. Debattista & SellwoodI l200(i lAthanassoulal l2003[ 


I Valenzuela & Klvpin I2OO3I) . The stellar disc could lose be¬ 
tween a few % and 40% of its angular momentum. Depending 
on the rate at which the stellar disc loses its angular momen¬ 
tum, the bar grows quite differently. The lack of dark halo in our 
simulations thus has mainly three consequences on our work: 

1) The evolution timescale could be different from simula¬ 
tions with a dark halo, leading us not to draw any conclusion 
about the evolution of ??, for instance. 

2) The shape and strength of a bar are also affected 
by th e presence of the live halo jAthanassoula & Misiriotisl 
I2OO2I) . The central co ncentration of the halo is thus a k ey pa¬ 


I il2nn2ll 


stress 


rameter. For instance. lAthanassoula & Misrrioti^ 
that more centrally concentrated halos lead to rectangular-like 
bars. However, the photometric and morphological properties 
(shape, ellipticity, Qt,, etc.; cf. also Paper I) of our simulated 
bars are very representative of real stellar bars. For instance, 
we are able to reproduce a rectangular-like bar (Run B, Fig.|^. 

3) The main results that could be affected by the absence 
of a dark halo in our simulations are those related to dynamical 
properties, i.e. the correlation of the position of resonances with 
bar-radii estimators (Sect. 16.H . This is a major concern since 
there is a debate on whether a dark halo increases 'R or not, 
leading to significantly small er bars compared to corotation. 
Debattista & Sellwood (120001) report high values of R for a 


numbe r of their simulations. However. IPebattista & SellwoodI 
1I2OOOI) estimate bar radii using in part the F 3 criterion that is, for 
our simulations, quite correlated to UHR. This criterion gives 
rather short bars, since it is mainly influenced by the part of the 
bar sustained by xi orbits. This might explain the high R va lues 
in part. For the simulations b v lO’Neill & Dubinskilil200^ and 


^^enzuela_&K 1 ^ 2 in i2003), R remains in the range 1.2-1.7. 


Valenzuela & Klvpinl ) I2003^ mainly use a criterion based on 
surface density (not tested in the present paper), which seems 
to giv e similar results than F 3 , whereas lO’Neill & Dubins^ 
1I2OO3I) use an average of three methods (including £3 and F 3 ). 

Our results should thus be conhrmed with simulations that 
include a live dark halo. However, the presence of a dissipa¬ 
tive component that can also exchange angular momentum and 
mass with the collisionless components through gravitational 
interaction, star formation, and feedback processes must have 
some influence on the global angular momentum exchanges. 
Therefore, realistic models should now include both dark mat¬ 
ter and a dissipative component. This will be published else¬ 
where. 


7 . Conclusions 

We have used N-body simulations, including stars, gas, and 
star formation, that were photometrically calibrated in B and 
K wavebands to compare the bar lengths determined by six ob¬ 
servational criteria (hve commonly applied to observations or 
simulations and a sixth introduced in this paper). These crite¬ 
ria are based on ellipse htting (Ei to E 3 ) and Fourier analysis 
(Fi to F 3 ) of the surface brightness (or mass surface density). 
The bar-length estimates were also compared to the location of 
resonances. 


We have obtained the following main results for each crite¬ 
rion in turn: 

- El dWozniak & Pierc^ Il99ll) . the radius of the maximum 
ellipticity, gives the shortest bar lengths. It clearly underes¬ 
timates the real bar radius. Moreover, it is not linked to any 
dynamical characteristics of a bar (e.g. resonance). 

- E 2 is measured at the end of the PA plateau. It can be used 
to approximately locate the UHR. 

- E3, th e radius of the minimum of ellipticity dWozniak et all 
1 19951) . gives the greatest bar lengths. In general, it is less 
sensitive to dust absorption than other criteria. This is the 
only criterion, amongst the six criteria tested, that is corre- 
lat ed with the corotat ion. 

- Fi dOhta et aklfT^ is not always dehned, since the bar- 
interbar brightness contrast (C) could remain under 2 even 
for strong bars. 

- F 2 dAguerri et al.Ll200(3) . also based on C, is the least noisy 
estimator. However it gives contradictory results for our 
simulations. 

- F 3 , widely used to determine the bar length in N-body 
simulations (e.g. IPebattista & Sellwoo(3l200r)l) . gives bar 
lengths comparable to E 2 results. Like E 2 , it gives quite a 
good approximation of the UHR location. 

In general, for projection angles i between 15 and 70°, the 
bar length is underestimated with respect to its face-on value, 
apart from using Ei. However, for i < 30°, errors remain below 
20%. Measurements also depend on the resolution but errors 
never exceed 15% when resolution is downgraded by a factor 5. 

We moreover conhrm that ellipticity prohles of B and K- 
band images can be very different. They also differ from ellip¬ 
ticity prohles of the mass density. This leads to a clear depen¬ 
dence of the maximum ellipticity on the waveband. However 
the radius of the maximum determined by Ei is less affected. 
Within a typical error of 10%, bar lengths determined with 
other criteria are not colour dependent, in the absence of dust. 
The situation is different when dust absorption is taken into ac¬ 
count. Criteria giving the smallest bar-length estimate are much 
more affected than others. This is the case of Ei in particular. 

For each simulation, the maximum ellipticity decreases 
with time, as do the bar strengths Q/,, but not at the same rate. 
The use of the maximum ellipticity to estimate the bar strength 
could thus lead to severe errors, especially in a comparison of 
one galaxy to another. 

Accurate determination of the bar length is crucial in sev¬ 
eral observational or theoretical analyses. For instance, the ratio 
R - Rcg/Rbar could increase from 1 (fast bar) to 1.4 (slow bar) 
just by using E 2 or F 3 instead of E 3 . However, being correlated 
to the corotation radius, E 3 cannot be used to dehne the end of a 
bar. Using E 2 or F 3 to dehne the end of the bar has physical im¬ 
plications since these criteria point out the UHR. Indeed, inside 
the UHR, the bar is mainly sustained by xi or other families 
of orbits elongated along the bar major-axis. However, there 
is clear evidence that what is commonly called a bar extends 
a bit outside the UHR because the region between UHR and 
CR is populated by higher resonant families of orbits, as well 
as semi-chaotic and chaotic orbits that contribute to the shape 
(sometime rectangular-like) of the bar. The dehnition of the bar 
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length (hence the criterion used) should thus depend on the ap¬ 
plication. 
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